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ABSTRACT 



Aims. We measure and study the evolution of the UV galaxy Luminosity Function (LF) at z=3-5 from the largest high-redshift survey 
to date, the Deep part of the CFHT Legacy Survey. We also give accurate estimates of the SFR density at these redshifts. 
Methods. We consider ~ 100000 Lyman-break galaxies at z ~ 3.1, 3.8 & 4.8 selected from very deep ugriz images of this data set 
and estimate their rest-frame 1600 A luminosity function. Due to the large survey volume, cosmic variance plays a negligible role. 
Furthermore, we measure the bright end of the LF with unprecedented statistical accuracy. Contamination fractions from stars and 
low-z galaxy interlopers are estimated from simulations. From these simulations the redshift distributions of the Lyman-break galaxies 
in the different samples are estimated, and those redshifts are used to choose bands and calculate k-corrections so that the LFs are 
compared at the same rest-frame wavelength. To correct for incompleteness, we study the detection rate of simulated galaxies injected 
to the images as a function of magnitude and redshift. We estimate the contribution of several systematic effects in the analysis to test 
the robustness of our results. 

Results. We find the bright end of the LF of our u-dropout sample to deviate significantly from a Schechter function. If we modify 
the function by a recently proposed magnification model, the fit improves. For the first time in an LBG sample, we can measure down 
to the density regime where magnification affects the shape of the observed LF because of the very bright and rare galaxies we are 
able to probe with this data set. We find an increase in the normalisation, <p*, of the LF by a factor of 2.5 between z ~ 5 and z ~ 3. 
The faint-end slope of the LF does not evolve significantly between z ~ 5 and z ~ 3. We do not find a significant evolution of the 
characteristic magnitude in the studied redshift interval, possibly because of insufficient knowledge of the source redshift distribution. 
The SFR density is found to increase by a factor of ~2 from z ~ 5 to z ~ 4. The evolution from z as 4 to z ~ 3 is less eminent. 

Key words, galaxies: high-redshift - galaxies: luminosity function - galaxies: evolution 



1. Introduction 

The formation and evolution of galaxies rank among the big 
questions in astronomy and still await a complete explanation. 
According to current theory, the formation of dark matter haloes 
by gravitational i nstabilities is an ess ential first step in the forma- 
tion of galaxies dEggen et alj|1962|) . Stars are be lieved to form 
when gas cools at the centres of these haloes (IWhite & Reesl 
and make up the part of the galaxy that we can observe. 
A number of physical processes strongly affect this baryonic 
mass assembly, like the hydrodynamics of the gas, feedback 
processes by supernovae and stellar winds, possibly magnetic 
fields, the role of AGN, or the effects of galaxy-galaxy inter- 
actions and mergers. For these reasons the modelling of galaxy 
formation depends on many free parameters and is not very well 
constrained. 



Over the past decade the high redshift universe has be- 
come accessible observationally through the use of photomet- 
ric techniques. By detection of the spectral discontinuity due 
to the redshifted Lyman-break in a multi-wavelength filter set, 
large and clean sa mples of high redshift s t ar-forming galax - 
ies can be selected dSteidel et al.ll 19961 1 1999t lGiavaliscoll2002h . 
with low amounts of contamination. These samples can be used 
to study several properties of the early universe. For example, 
by measuring the correlation function of these Lyman-break 
Galaxies (LBGs) and comparing it with the correlation of dark 
matter, the characteristic mass of their haloes can be deter- 



mined (e.g. iGiavalisco & Dickinson! 1200 it lOuchi et all l2004bl 
iHildebrandt et alj|2005l I2007L 120091) . Hubble Space Telescope 
observations of LBGs are use d to study how cert ain morpho- 
logical types evolve with time dPirzkal et al.ll2005l). A study of 
the evolution of the UV Luminosity Function (LF) ([Steidel et al.l 
119961 119991: ISawicki & Thompsonll2006l iBouwens et al.ll2007l) . 



* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de l'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is based in part on data products 
produced at TERAPIX and the Canadian Astronomy Data Centre as 
part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- 
rative project of NRC and CNRS. 



which is the measure of the number of galaxies per unit vol- 
ume as a function of luminosity, is another fundamental probe in 
galaxy formation and evolution, because of its close relation to 
star formation processes. 

Several techniques can be used to estimate the star formation 
rate (SFR) in galaxies, mostly based on the existence of massive, 
young stars, indicative of recent star formation. A commonly 
used way to probe the existence of massive stars is the Ha- lu- 
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minosity dKennicut til 19831) . because Ha photons originate from 
the gas ionized by the radiation of massive stars. A second star 
formation indicator is the infra red (IR) luminosity originating 
from dust continuum emission (lKennicutjll998t iHirashita et al.l 
120031) . The absorption cross section of dust is strongly peaked in 
the UV, and therefore the existence of UV emitting, i.e. mas- 
sive, stars is probed indirectly. Thirdly, the UV continuum is 
used as a star formation probe, with the main advantages being 
that the UV-emission of the young stellar population is directly 
observed, unlike in H-a and IR studies. Furthermore, this tech- 
nique can be applied from the ground to star-forming galaxies 
over a wide range of redshifts. Hence, it is sti ll the most powerfu l 
probe of cosmological evolution of the SFR (Ma dau et al.l 19961) . 
However, information about the initial mass function (IMF), and 
particularly the extinction by dust are required to estimate the 
total star formation rate. 

In this paper we estimate the UV LF at redshifts z=3- 
5 from the Canada-France-Hawaii-Telescope Legacy Survey 
(CFHTLS) Deep, a survey covering 4 square degrees in four 
independent fields spread across the sky. Since our samples, at 
different redshifts, are all extracted from the same dataset, this 
gives an excellent opportunity to study a possible evolution of 
the LF in this redshift interval. Several systematic effects that 
need to be considered when comparing results at different red- 
shifts from different surveys (e.g. source extraction, masking, 
PSF-modelling, etc.) can be avoided when the different redshift 
samples are extracted from the same survey. Due to the large vol- 
umes we probe with our 4 square degree survey, the influence of 
cosmic variance on the s hape of the estimated LF is negligible 
dTrenti & Stiave lli 2008), as we expect cosmi c variance to af- 
fect o ur number counts only at the 1-5% level (ISomerville et al.l 
12004 . We can study the bright end of the LF with unprecedented 
accuracy, as these objects are rare and we are able to measure 
down to very low densities. This allows us to study the effect that 
magnifi cation has on the observed distribution (see recent re- 
sults by I Jain & Limall20Tol) . and study a possible deviation from 
the commonly used Schechter function. Furthermore, given the 
depth of the stacked images, we can probe the faint end of the 
luminosity function with comparable pr ecision as the deepest 
groun d based surveys have done before (ISawicki & Thompson! 
120061) . 

The structure of this paper is as follows: In Sect. [2] we de- 
scribe the data set we use, the LBG selection criteria as well as 
the simulations that lead to the redshift distributions and contam- 
ination fractions. In Sect.[3]the survey's completeness and the ef- 
fective survey volumes are estimated. In Sect.|4]we proceed with 
the resulting estimated LFs, present the best-fitting Schechter pa- 
rameters, and show how a simple magnification model can sig- 
nificantly improve the quality of the fit. The UV luminosity den- 
sities (UVLD) and SFR densities (SFRD) are estimated based 
on the measured LFs. We also elaborate on the robustness of our 
results. In Sect. [5] we compare these to previous determinations 
of the UV LF and SFRD from the literature. In Sect.[6]we finish 
with a summary and present our conclusions; 

We use the AB magnitudes system (lOke & Gunnl 1 1983b 
throughout and adopt ACDM cosmology with £2 m = 0.3, 
f2 A =0.7 and H = 70kms-' Mpc" 1 . 



2. Data & Samples 

2.1. The CFHT Legacy Survey Deep 

For this work we make use of publicly available data from the 
CFHTLS Deep, a survey using MegaCam mounted at the prime 



focus of the CFHT which covers four independent fields of 1 
square degree each. Images are taken in t he filters ugriz and are 
pre-p rocessed using the Elixir system dMagnier & Cuillandrd 
12004. Image reduction is done using t he THELI pipeline 
(lErben et al.l 120051 iHildebrandt et al J 12003) . leading to approx- 
imate 5cr point source limits of 27.5, 27.9, 27.9, 27.7 and 26.5 
for the ugriz bands, respectively. The limits for each of the fields 
lie within 0.2 mag fro m these average values. 

Source Extractor dBertin & Arnoutsll 19961) is used to create 
a multi-colour catalogu e. Total ;-ban d magnitudes are measured 
in Kron-like apertures (iRronl IT980h using SExtractor's AUTO 
magnitudes. Every image is smoothed by convolution with a 
Gaussian filter to match the seeing of the image with the worst 
seeing value (typically the w-band). These corrections are typi- 
cally small - all bands have a seeing below 1 arcsec. The dual- 
image mode of SExtractor is then used with the unconvolved 
/-band for source detection and isophotal magnitudes from the 
convolved bands to estimate colours. An extensive d escription of 
the da ta re duction and catalogue cre ation is given in lErben et al.l 
(120091) and lHildebrandt et al.l (l2009h . 



2.2. u-, g-, and r -dropout samples 

Clean samples of u-, g-, and r-dropouts are selected based on the 
following selection criteria dHildebrandt et al. 20091): 

M-dropouts : 1.0 < (u - g) A -1.0 < (g - r) < 1.2 A 

1.5 -(g-r)<(u-g)- 0.75, (1) 

g-dropouts : 1.0 < (g - r) A -1.0 < (r - i) < 1.0 A 

1.5 -(r-i)<(g-r) -0.80, (2) 

r-dropouts : 1.2 < (r — i) A —1.0 < (i — z) < 0.7 A 

1.5 • (i-z) <(r- 0-1.00. (3) 

Furthermore, it is required that all LBGs have a S Extractor 
CLASS.STAR parameter of CLASSJSTAR < 0.9, that g- 
dropouts are not detected in u, and that r-dropouts are neither 
detected in u nor in g. Note that the colour sele ction criteria of 
the M-dropout sample in lHildebrandt et all (l2009h are pretty con- 
servative. We relaxed the u — g cut slightly to m a ke the selec- 
tion m ore comparable to e.g. Steidel et al.l d2003l) . ISteidel et al.l 
d2004 and lSawicki & Thompson! d2006i 

This results in the selection of 50880 u-, 36226 g-, and 11411 
r-dropouts in total over the four fields. Their magnitude distri- 
butions are shown in Fig.[T] Note the differences in depth of the 
individual fields. 



2.3. Redshift distributions & Contamination fractions 

The majority of the selected sources is too faint to make a spec- 
troscopic redshift determination possible, and the brighter can- 
didates havejio^_sjje£^sco£ically been observed yet. For this 
reason Hildebrandt et al. I <l20Q9t) estimated the redshift distribu- 
tions by means of photometric redshifts and simulations. 

Throughout this paper we will use the mean redshift 
values from simulat i ons b ased on synthetic templates by 
iBruzual A. & Chariot! (1 19931) . being < z >= 3.1, 3.8 and 4.7 for 
the m-, g-, and r-dropouts respectively. We estimate the uncer- 
tainty in the mean redshifts to 0.1 for the three dropout samples. 

In order t o address the amount o f contamination in our 
LBG samples, Hil debrandt et al.l (120091) consider the possibili- 
ties of stars and low-z galaxies scattering in the selection boxes. 
Galaxies are simulated base d on templates from the library of 
IBruzual A. & Charloll (1 19931) . while the colours of stars in the 
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Number counts selected u— , g— and 7 — dropouts from the CFHTLS Deep 
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Fig. 1. From left to right the number counts of the u-, g-, and r-dropouts in the CFHTLS Deep, as selected by the colour criteria of 
Eqs.HU 
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Fig. 2. Contamination fractions of stars and low-z galaxies in the dropout samples. Blue ^-symbols connected by a dashed line show 
the stellar co ntamination fraction base d on a galactic model. Black squares connected by the solid line show the total contamination 
fraction from lHildebrandt et all (120091) . 



fields are estimated based on the TR1LEGAL galactic model 
dGirardi et al.ll2.0Q5h . Contamination fractions are shown graphi- 
cally in Fig. [2] 

Stellar contamination is negligible for the brighter objects, 
as the selection boxes steer away from the stellar locus. For faint 
objects in the u-, and g-dropout samples, the stellar contamina- 
tion increases as a result of photometric scatter. 

The contamination by low-z galaxies is negligible in the u- 
dropout sample, as the Lyman break for a z ~ 3 galaxy is still 
blue-ward of the z ~ Balmer/4000A break. For higher red- 
shifts this ceases to be the case, so that the g-, and r-dropout 
samples suffer from a significant contamination fraction at the 
bright end, where the LBG population is sparse. Faint low-z ob- 
jects are likely to scatter into the selection box in each of the 
samples, so that the contamination fractions are increased here. 

Since the contamination fractions are low at the bright end 
of the w-dropouts, we have the potential to probe the LF to very 
low LBG galaxy densities. We inspect the 80 brightest objects 
(>ab < 23.2) in the M-dropout sample by eye and remove obvious 
spurious sources, 30 in total. A strong, but certainly slightly sub- 
jective, rejection criterion is the size of the region in which the 
flux is measured, i.e. the half-light radius (13 objects rejected). 
Sources that are clearly blended by a bright neighbour (2 ob- 
jects), as well as sources that have a too large apparent size (3 
objects). Also an asteroid track has been removed. 



In the g-, and r-dropout sample we do not probe these low (< 
0.01 mag -1 arcrrrkr 2 ) LBG galaxy densities, since these points 
are unreliable due to the high contamination fractions from low- 
z objects (see Fig. |2j. This prohibits a study of the bright end of 
the LF for the g-, and r-dropouts until the nature of the individual 
sources has been verified. 



3. Analysis - Survey completeness 

We use a detailed modelling approach to estimate the complete- 
ness of the survey as a function of magnitude and redshift, for 
each of the dropout samples. We add artificial objects to our im- 
ages, with colours representative of star-forming galaxies, and 
try to recover them following the same source extraction and 
colour selection criteria as for the real data. We investigate how 
the increasing scatter in the colours for fainter objects influences 
the completeness as a function of magnitude. Furthermore, one 
expects that, for fainter dropout objects, the redshift distribution 
of these objects broadens due to the same effect. Hence we will 
also model this as a function of redshift. 

We describe our fiducial model SEDs, sizes and shapes of the 
simulated objects below. The assumptions we make, are tested in 
Sect.l4.3lto estimate the robustness of the results. 
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3.1. Model galaxies 



The iBruzual A. & Charlotl (I1993I) stellar population synthesis 
library is used to set up our fiducial galaxy SED model; a 
100 Myr old g a laxy t emplate with constant star formation. A 
iMiller & Scalol d 19791) IMF is assumed. The optical depth of 
neutr al hydrogen, as a function of redshift, is modelled accord- 
ing to lMadaulIT995l) . 

This template i s redd ened by the starburst extinction law 
from lCalzetti et al.l d2000l) . with a distribution in E(B - V). This 
distribution we choose such that the UV-continuum slopes that 
we measure from the data are matched when we use our fiducial 
template as a base. To measure the UV-continuum slopes, we 
use a colour redward of 1600A rest-frame, i.e. the r - i colour 
for the M-dropouts and the i—z colour for the g-dropouts. For the 
r-dropouts we can not perform a similar measurement because 
we do not have observations in a band redward of the z-band. 
Therefore we will use the same distribution of dust as we find 
for the g-dropouts. 

For the M-dropouts we find that a uniform dust distribution 
with reddening between with 0. 1 < E(B - V) <0.4 gives a good 
fit to the data, see Fig. [3] For the g-dropouts we measure a 
larger spread in UV-continuum slopes, and find a reasonable fit 
when using a uniform dust distribution with reddening between 
0.0< E(B - V) <0.5. It should be noted that the age and amount 
of dust attenuation of the template are highly degenerate, so that 
different combinations of these parameters fit the UV-continuum 
slopes in the data. It is especially important to correctly match 
the distribution of the UV-continuum slopes from the data with 
the model galaxies, since an increase (decrease) in the age of 
the galaxy model template has a similar effect on the LF as an 
increase (decrease) of the amount of dust. We will elaborate on 
this in Sect. 14.31 If we measure the UV-continuum slopes in our 
data for different magnitude bins, we do not find a significant 
evolution. Therefore we will use a distribution of dust attenua- 
tion in our simulations that is independent of magnitude. 

LBGs have typica l half-light radii of r\/2 ~ 0.1"-0.3" 
(iGiavalisco et al.Nl 9961) and thus are unresolved by the CFHT 
and can be treated as point sources. As the size of the PSF in the 
CFHT images is strongly position dependent, we have to adapt 
the injected sources accordingly. We parametrize the PSF by a 
Moffat profile, 



/ = /» 



1 + (2 



m. 



(4) 



in which B and Rq are the parameters that we adapt to adjust the 
shape and size of the profile respectively. Iq represents the flux 
normalisation. 

In order to measure the PSF as a function of position for the 
different filters and fields, we first select several hundred stars 
based on their magnitudes and half-light radii, in each field and 
each filter. We measure the 50% and 90% flux radii of these 
stars, rso and rgo, using SExtractor. The ratio of these flux radii 
uniquely determines the Moffat-/? parameter, which, in combi- 
nation with either one of the flux radii, gives the Moffat profile 
radius Rq for each star. We find that the Moffat-/? parameter is 
fairly constant over the fields and filters {B a 4.0) and only Rq 
changes significantly. To model Rq as a function of position, we 
fit a 2-dimensional polynomial function to the SExtracted values 
for T5o. This constrains the PSF size on every position, for every 
field and filter. 

We find that there is a ~30% difference in Rq between the 
image center and boundaries. As this, in our ground-based wide- 



CFHTLS field 1 
10Myr cont, E(B-V)=0.15 
100Myr cont, E(B-V)=0.15 . 
200Myr cont, E(B-V)=0.15 
100Myr cont, dust distribution. 




Fig. 3. The distribution of UV-continuum slopes for the u- 
dropouts as measured by the r - i colour from the CFHT data 
(black) compared with the outcomes from simulations. The 
colour is measured in a part of the spectrum without strong fea- 
tures, redward of 1600A. Both an older template and a higher 
amount of dust result in a redder UV colour. The amount of dust 
we need to add therefore depends on the age of the base tem- 
plate. For the M-dropouts we find that a uniform distribution of 
dust with 0. 1 < E(B - V) <0.4 gives a good fit to the data, when 
the template is 100 Myr old with a constant amount of star for- 
mation (purple). 



field survey, is by far the dominating effect in the apparent sur- 
face brightnesses of our sources, we do not assume an intrinsic 
size distribution for the sources in the main simulations of this 
paper. 

3.2. Eddington bias 

If we consider a certain intrinsic magnitude distribution of galax- 
ies, the recovered distribution after source extraction and colour 
selection will look different due to statistical fluctuations in the 
measurement. In our anal ysis we attempt to correct for this ef- 
fect called Eddington bias dEddingtonl 1 9 1 3t lTeerikorpil2004l) . It 
is hard to estimate this bias analytically because the size of the 
magnitude scatter is an increasing function of magnitude. Also, 
when approaching the completeness limit of the survey, only the 
brighter objects will be detected. What generally happens is that 
intrinsically faint objects will on average look brighter than they 
are. 

The Eddington bias can be estimated by comparing the in- 
trinsic magnitudes of the injected sources to the recovered mag- 
nitudes. Such a comparison is shown in Fig. [4] stressing the im- 
portance of that effect for faint magnitudes. 

We want to inject sources with the same magnitude distribu- 
tion as the intrinsic distribution, to get an unbiased result. As the 
bias is expected to be largest for fainter sources, it is especially 
important to correctly model the faint-end of the intrinsic distri- 
bution. Here we adopt a LF that is consistent w ith the deepest 
LBG survey, conducted bv lBouwens et al ] (120071) . 



3.3. Source injection and recovery 



Following this adopted intrinsic distribution we inject 20000 
sources in each of the images, for 60 equal redshift steps be- 
tween z = 0.0 and z = 6.0. To verify that the injected sources do 
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Fig. 4. The difference between recovered and injected source 
magnitudes as a function of recovered magnitude for the r- 
dropouts. Similar trends appear in the other dropout samples. 
The error bars reflect the scatter from the four fields. A sim- 
ilar effect as shown ar ound m lec = 25 . 3, nam ely an increase of 
m lec - rrtinj , is found bv lBouwens et alJ (120061 Fig. 1 9), where ob- 
servations were used rather than simulations. 



not significantly influence each other by blending, nor that the 
background is influenced significantly, we perform the following 
tests. We inject the same 20 000 sources in 4 stages, 5000 sources 
each, and do a third analysis where we inject 100000 sources in 
total. In Fig.[5]the recovered fractions of sources that also satisfy 
our g-dropout criteria are shown as a function of magnitude, for 
one particular redshift step. Only for faint magnitudes does the 
100000 curve deviate from the other ones, which are identical 
in this regime. In Fig.|6]the distribution of recovered g-dropouts 
with an intrinsic magnitude of m = 25.0 is shown as a func- 
tion of recovered magnitude. We conclude that the injection of 
20 000 sources does not influence the images such that the pho- 
tometry would be perturbed significantly. A similar behaviour is 
expected for the u- and r-dropout samples. 

The clustering of LBGs is not taken into account. We as- 
sume this effect to be insignificant for estimating complete- 
ness, as the correlation l ength, which is typically around 5 Mpc 
dHildebrandt et al.l l2009). is very small compared to the survey 
volume. Therefore we spread our simulated sources uniformly 
over the images. 
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Fig. 5. The recovered fraction of injected sources that also satisfy 
the g-dropout criteria, as a function of /-band magnitude, for 1 
redshift step. Three curves are given for different amounts of 
injected sources, to see whether, and how, the presence of these 
sources influences the photometry of the analysis. 
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Fig. 6. The distribution of measured source magnitudes for a 
population of injected simulated sources with intrinsic /-band 
magnitude=25.0, for 1 redshift step. Three curves are given for 
different amounts of injected sources, to see whether, and how, 
the presence of these sources influences the photometry of the 
analysis. 



3.4. Effective volumes 

Next we define the function p(m, z) to be the number of sources 
recovered with an observed magnitude in the interval [m; m + 
Am], and are selected as dropouts, divided by the number of in- 
jected sources with an intrinsic magnitude in the same interval 
[m;m + Am] and a redshift in the interval [z;z + Az]. Note that 
the definitio n of p(m, z) is slightly differen t compared to the one 
used in e.g. ISawicki & Thompson! d2006l) . as they do not take 
Eddington bias into account. In our definition, p(m, z) could po- 
tentially be > 1 as a result of this bias correction. 

The effective volumes ( V e ff) of our survey are given by 

CdVc 

V eS (m)=Af ——p(m,z)dz, (5) 



where A j is the field area in square arcminutes, and -p is the 
comoving volume per square arcminute, which depends on the 
adopted cosmology. 

The magnitude is measured in the r-, /-, and z-bands for the 
u-, g-, and r-dropout samples, respectively. These bands probe 
flux at approximately 1600 A rest-frame of the sources at the ex- 
pected mean redshifts, so that only a minor k-correction will be 
sufficient to compare the results for the different epochs directly, 
see the upper panel of Fig. [7] We transform the apparent mag- 
nitudes to absolute magnitudes and perform a k-correction to a 
rest-frame wavelength of 1600 A using the mean redshifts for 
each of the dropout samples, i.e. we assume all sources to be at 
the same redshift. 
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100 Myr, E(B-V)-Q.25 




Table 1. The estimated LFs from the CFHTLS Deep. 



Fig. 7. Upper panel: k-correction to 1600 A for the MegaCam 
riz filters, as a function of redshift based on the 100 Myr old 
continuously star-forming template with a dust attenuation of 
E(B - V)=0.25. The average redshifts from the simulations and 
their uncertainties are represented by horizontal error bars. This 
leads to a corresponding error in the k-correction. Lower panel: 
Shifts from apparent magnitudes in the MegaCam riz filters, to 
absolute magnitudes at 1600A. The Distance Modulus and k- 
correction are taken into account. The uncertainties in the aver- 
age redshifts of the samples lead to uncertainties in the absolute 
magnitudes. 



The uncertainties in the mean redshifts[3 are expected to be 
approximately 0.1 for the three dropout samples, as argued in 
Sect. 12.31 This leads to uncertainties in both distance modulus 
and k-correction, resulting in a systematic error in the absolute 
magnitudes of our estimated LF, see the lower panel of Fig. [7] 
The final uncertainties are about 0.07, 0.05 and 0.04 in the ab- 
solute magnitude for the u-, g-, and r-dropout samples, respec- 
tively. 

4. Results 

4.1. The UV Luminosity Function at z=3-5 

After dividing the number counts by V e s, which corrects for in- 
completeness and Eddington bias, and subtracting the distance 
modulus and the k-correction from the apparent magnitudes, we 
obtain the LF in absolute magnitudes at 1600A. 

Results from the four fields are binned to Amag=0.3, com- 
bined, and shown in Fig.[8]and TableQ] Uncertainties in the mag- 
nitude direction are due to uncertainties in the redshift distribu- 
tion of source galaxies. The four independent analyses of the 
fields allow us to estimate field-to-field variations for each of the 



1 Note that we make use of two different redshift distributions in our 
analysis. To estimate both the k-correction and the effective volumes we 
use the distribution given by the simulations describe d in Sect. [3] while 
we us e the distribution from the simulations from Hildebrand t et alj 
d2009l) to shift the LF in the magnitude direction. The redshift distri- 
bution from the latter simulations are expected to be more reliable be- 
cause they take a wide variety of template galaxy models into account, 
but we cannot use them for effective volumes & k-corrections because 
of computational constraints. Rather we have to simulate those with a 
single template. 
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data points. Vertical error bars reflect either this uncertainty, or 
the Poisson noise term, whichever is largest. Usually the field- 
to-field variance dominates. As a consequence of the way these 
are computed, Poisson noise is a lways taken into account. 

We fit a Schechter function (ISchechterlll976h to the binned 
data points, 



<p(M)dM = 0.4 ln( 10)0* 10 



0.4(ff+l)(M*-M) 



exp(-10' 



OA(M'-M) 



), (6) 



with M* being the characteristic magnitude, a being the faint- 
end slope, and (p* being the overall normalisation. 

Using x 1 statistics on a three dimensional grid of 500 3 dif- 
ferent Schechter parameter combinations, we find the minimal 
value 0t^ in ) for each of the dropout samples yielding the best fit 
values. To estimate the errors in the fitted parameters, we project 
the 3-dimensional distribution of x 2 to 3 planes by taking the 
minimum^ 2 along the projected dimension. In Fig.|9]the 68.3% 
and 95.4% confidence levels are shown, which correspond to a 
A^ 2 =2.3 and 6.17 with respect to^f^. In Table [2] we give the 
68.3% confidence levels on each individual parameter, corre- 
sponding to Ay 2 =1.0. 

Note, however, that this error estimate assumes Gaussian 
errors, and that the errors on the data points are independent. 
Especially for the w-dropouts this is probably not the case. The 
normalisation of the LF seems to be systematically slightly dif- 
ferent for each of the fields (see also Fig.[T|), giving a systematic 
uncertainty in 0* of about 30%. For the M-dropouts this effect 
is largest because of a slightly more uncertain flux calibration 
in the M-band compared to the g- and r-bands. The effective fil- 
ter throughput changes with time in the UV as the atmosphere is 
changing, and also the camera is less sensitive in this wavelength 
regime resulting in larger shot-noise. 

4.1.1. Magnification contribution at low densities 

Due to inhomogeneities in the matter distribution between dis- 
tant sources and the observer, paths of photons get slightly per- 
turbed. This results in a distortion of the shape and a magnifi- 
cation of distant sources. When a source is magnified by a fac- 
tor fi, the flux gets boosted by the same amount. One can re- 
late an intrinsic luminosity distribution to an observed distribu- 
tion if the magnific a tion distribution is k nown, as was shown by 
Ijain & Lima! ( 2010h . lHiibert et alJ d2007l) estimate magnification 
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Fig. 8. The LFs of LBGs in the CFHTLS-Deep fields. Data points and best-fitting Schechter functions are shown for the w(blue)-, 
g(green)-, and r(red)-dropouts. For legibility we applied a small offset on the x-axis values of the g-, and r- dropouts. The dashed 
blue curve shows the best fitting Schechter function for the M-dropouts after magnification effects have been included, as described 
in Sect. l4~TTT1 
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Fig. 9. The 68% and 95% likelihood contours for different 
Schechter parameter combinations. Shown are the results for the 
M-dropouts at z ~ 3 (blue, solid), the ^-dropouts at z ~ 4 (green, 
dots), and the r-dropouts at z ~ 5 (red, dash-dots). The u-dropout 
contours represent the best Schechter parameters when we in- 
clude magnification effects, as described in Sect. I4.1.1I 



distributions for different source redshifts by shooting random 
rays through a series of lens planes created from the Millennium 
Simulation. The width of the magnification distribution is found 



Table 2. A comparison between the best fitting Schechter pa- 
rameters and their 68% confidence intervals for the u-, g-, and 
r-dropouts. 



Sample 


M* a 
lvl uv 


P 

[lO^Mpc" 3 ] 


a 


^ 2 /dof 


u 

g 
r 

u b 


-20.94^-^ 
-20.84+™ 
-20.94+°]? 
-20.84*5 


1 70+0.5 1 
1 36+°- 23 

1. j -0.20 

0.83!;^ 

9 1 1 +0.63 


-1 65+ ,u - 
-1 56+°- 08 
-1 65+ - 09 

1,UJ -0.08 

-1 60+°- 14 


0.52 
0.36 
0.19 
0.41 



Due to uncertainties in the redshift distributions there is an addi- 
tional error component of 0.07, 0.05, and 0.04 in the estimated M^ v 
for the 11-, g-, and r-dropout samples, respectively. 
Best-fitted Schechter parameters for a model where the function is 
modified with a magnification distribution. 



to increase with increasing source redshift, and the peak posi- 
tion of the distribution decreases slightly with increasing source 
redshift. 

Magnification can account for a strong deviation from a 
Schechter function where the slope of the i n trinsic luminosity 
distribution is very steep, see iJain & Lima! (l2010l) . We stress 
again that we measure the LF from a volume that is much 
larger than has been used before. The bright end of our g-, 
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and r-dropout samples suffers from increasing amounts of con- 
tamination. Only for the M-dropouts we probe the distribu- 
tion of M-dropouts at the bright end down to a density of 
10 mag 1 Mpc~ 3 . We use the magnification distribution for a 
source redshift of z = 3.1 that was kindly provided by Stefan 
Hilbert to improve our model to the data. 

Writing the LF as a function of magnitude, we use the 
following equation to correct the Schechter function for mag- 
nification effects. I t is equivalent to the expression used by 
|Jain& Lima! (120101) . 



4>( m obs) 



■I 



d/j, PQj.) <p*(m obs + 2.5 logOu)), 



(7) 



where <f>* is the Schechter function defined by Eq. [6] The new 
function yields a slightly better fit to the bright end of the LF, 
reducing the formal ^ 2 /dof from 0.52 to 0.41. Find the new 
Schechter parameters, together with their 68.3% confidence lev- 
els in Table |2 and their 2-dimensional 68.3% and 95.4% confi- 
dence contours plotted in Fig. [9] The best fitting function is the 
dashed curve in Fig. [8] 

As the bright selected M-dropouts are likely to be signifi- 
cantly magnified, we expect them to appear close to a massive 
foreground galaxy or group of galaxies that acts as a lens. We 
inspected the brightest (r^s <23.2) M-dropouts by eye and find 
that this is indeed the case for many of them. Note that the 
model is still uncertain as the Millennium simul ation does not 
includ e baryonic matter, and assumes cr% = 0.9 (ISpringel et al.l 
12001. where recent estimates indicate a lower value around 0.8 
( Komatsu et al.l20l"oT> . Also, in the magnification probability dis- 
tribution the possibility of multiply imaged systems is ignored. 
To rule out contamination in the LBG sample at the bright end of 
the LF, the nature of each bright object has to be verified spec- 
troscopically. 

A statistically much better sample can be selected from the 
CFHT Legacy Survey Wide, consisting of 170 square degree 
imaging in ugriz of shallower depth. We leave this for future 
studies. 



4.2. The UV Luminosity Density and SFR Density 

Next we estimate the UV luminosity density (UVLD) at the dif- 
ferent epochs. To be able to compare our results with the results 
from previous studies, we will integrate the data points down to 
L = 0.3L! =3 , where L* =3 (= 9.4 x 10 28 erg s" 1 Hz" 1 at 1600A) is 
the characteristic luminosity of our M-dropout sample. Results 
are shown in the odd-numbered rows of Col. 3 in Table [3] 
However, for steep faint-end slopes of a < -1.6, more than 50% 
of the UV luminosity is expected to be emitted by lower lumi- 
nosity sources. Therefore we will make a second estimate of the 
UV luminosity density by integrating the best-fitting Schechter 
function over all luminosities. This results in 



puv = <p*L*T(a + 2), 



(8) 



where T is Euler's Gamma function. Although this full inte- 
gral of the LF depends strongly on uncertainties in the faint-end 
slope, we use it to provide an upper limit to the UVLD. The re- 
sults are shown in the even-numbered rows of Col. 3 in Table [3] 
The effective extinction in the UV is a strong function of the 
amount of dust. At these high redshifts (z > 3) the only esti- 
mate for the amount of dust is based on a measure of the UV 
continuum slope. Note however that there is a strong degeneracy 
between the age and the amount of dust in the template i f the 
rest-frame IR is not covered, see e.g. lPapovich et al1(l200ll) . 



Table 3. The puv and psfr for the different dropout samples. 
The first lines for each sample correspond to sums over the data 
points down to L — 0.3L*_ 3 while the second lines correspond to 
integrals over the best-fit Schechter functions. 



Sample 


Integral limit 


Puv [10 ift 
ergs -1 Hz 1 Mpc~ 3 ] a 


Psfr 

[M yr' Mpc- 3 ] ab 


u 


L > 0.3L* =3 
Schechter 


1.73 ±0.09 
4.41 


n i 0Q+U.U64 
ua, "-0.036 
0.154 


8 


L > 0.3L* =3 
Schechter 


1.07 ±0.03 
2.62 


0.07832® 
0.092 


r 


l > o.3l: =3 

Schechter 


0.80 ± 0.03 
2.19 


027 +0 013 

"•"^ ' -0.007 

0.038 



Due to uncertainties in the redshift distributions, there is an addi- 
tional error component of ~7%, ~5%, and ~4% in the estimated 
Puv and p S FR values for the u-, g-, and r-dropout samples, respec- 
tively. 

Corrected for dust e xtinction using the lum inosity dependent cor- 
rection factors from iBou wens et all d2009h . Systematic errors as a 
result from the age-dust degeneracy are also included. 



iBouwens et al.l d2009l) recently measured the UV-continuum 
slope of LBGs at high redshifts from deep HST data, from which 
the amount of dust obscuration could be estimated as a func- 
tion of LBG magnitude. The values they find at the characteristic 
magnitudes of our samples are E(B - V) — 0.15 forz - 3,4 , and 
E(B-V) = O.lOforthez ~ 5 sample. lBouwens et alj d2009t) find 
a decreasing amount of dust for fainter magnitudes. For consis- 
tency we will use the relationships they estimate to correct for 
dus t extinction in our dat a, and not the values from Sect. 13.11 

iMeurer et al. I (fl99l find a relation between the UV- 
continuum slope and the extinction by dust. IBouwens et al.l 
(120091) use this relation and find, upon integrating down to 
L = 0.3L* =3 , density correction factors of 6.0+} ^, 5.8+°*^, 
and 2.7^q ^ ( ) 7 for the three redshift samples, res pectively. Both 
rando m errors and systematic errors are quoted ( IBouwens et al.l 
I2009I 

We now convert the UV luminosity density into the star 
formation rate den sity, psfr, at the different epochs using 
dMadauet al.ll 1998b . 



L^S.OxlO 27 l-^-] ergs 1 Hz 1 . 
\Moyr- 1 ; 



(9) 



This relation assumes a 0. 1-125M0 Salpeter IMF and a constant 
star formation rate of > 100 Myr. The resulting estimates of psfr 
are shown in Col. 4 of Table [3] where we have corrected for dust 
extinction. In Sect. 15.41 we compare these estimates to values 
reported in previous studies. 

Note, however, that some sources like AGN, which might be 
included in our dropout samples, add to the total UV luminosity 
density in the Universe, though do not contribute to the SFRD. 

4.3. Robustness of our results 

Our fiducial template model is a 100 Myr old continuously 
star-forming galaxy with a uniform distribution of dust centred 
around E(B - V)=0.25. This dust distribution was chosen such 
that the distribution of UV-continuum slopes of the recovered 
simulated sources matches the distribution of UV-continuum 
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slopes in the real data (see Fig. [3}. We test some of the assump- 
tions we made in Sect. 13.11 by checking their influence on the 
final LFs. 

As a reference SED we use a 100 Myr old galaxy model 
with constant star formation and a single dust attenuation value 
of E{B - V)=0.25. We consider redder (bluer) templates by either 
increasing (decreasing) the age of the star-forming period, or in- 
creasing (decreasing) the amount of dust. In Fig.[10]the colours 
of these alternative templates, as they would be measured by the 
MegaCam ugriz filter set, are shown as a function of redshift. As 
the quality of the Schechter fit is high in all cases (^ 2 /dof < 1 .0), 
we present the differences by comparing the Schechter parame- 
ters: 

- The faint-end slope a depends on the colour of the spectral 
template. In the g-, and r-dropout samples, redder templates 
tend to give steeper a than bluer templates. The reason for 
this is as follows, with the g-dropouts as an example. Sources 
in the selection box have red observed g - i colours. For faint 
magnitudes (note that the reference magnitude is measured 
in the /-band), the g-band magnitude exceeds the limiting 
magnitude of the survey. Only a lower limit on g - r can then 
be given and the source moves out of the dropout selection 
box. The magnitude at which this happens depends on the 
g - i colour of the template. As the average g - i colour in 
the selection box is redder for red templates (Fig. ITOb . the 
detection rate of red sources is suppressed at faint magni- 
tudes. This argument holds for any red template in the g-, 
and r-dropout samples. Fig. [TO] indicates that the opposite 
effect happens in the M-dropout sample as the u — r colour 
is generally bluer (redder) for redder (bluer) templates in 
the selection box. This effect is indeed also inferred from 
the simulations. There is an additional effect due to the re- 
quirement that g-dropouts are not detected in u, and that r- 
dropouts are neither detected in g nor in u. This suppresses 
the detection of bright g-, and r-dropouts, especially at low 
redshifts. The bluer the u-i colour for the ^-dropouts (i.e. the 
bluer the template), the stronger this effect is. A similar argu- 
ment holds for the r-dropouts. Therefore the V e ff's are higher 
in the faint magnitude regime, so that the LFs are lower at 
this end. The ranges of best-fitted a values, for the template 
spectra we considered, are -1.82 < a < -1.38 for the u-, 
-1.82 < a < -1.42 for the g-, and -2.14 < a < -1.40 for 
the r-dropout sample. 

- The characteristic magnitude, M*, does not sensitively de- 
pend on the template spectrum chosen. The ranges of best- 
fitted M* values, for the template spectra we considered, are 
-20.88 < M* < -20.72 for the u-, -21.06 < M* < -20.74 
for the g-, and -21.26 < M* < -20.88 for the r-dropout 
sample. However, another systematic uncertainty in this pa- 
rameter is due to the unknown redshift distribution of source 
galaxie s, see Fig. Ill which depen ds on the mix of templates 
used bv lHildebrandt et alj J20091) . 

- The normalisation 0* decreases (increases) when the 
faint-end slope becomes steeper (shallower). The best-fit 
Schechter parameters move then in the direction of the de- 
generacy of the ellipse in the upper left part of Fig. [9] The 
ranges of best-fitted <f>* values [10 -3 Mpc -3 ], for the tem- 
plate spectra we considered, are 1.36 < <p* < 2.89 for the 
u-, 0.76 < (p* < 1.74 for the g-, and 0.39 <</>*< 1.12 for the 
r-dropout sample. 

Some studies (e.g. ISawicki & Thompsonl 12006) make use 
of a starburst template instead of a continuously star-forming 



model. The stellar population in a starburst template is older 
on average, and therefore the colours will be redder. However, 
for a template age of 100 Myr the difference in colours is very 
small. We compare the Schechter parameters that we measure af- 
ter using our reference model (i.e. a 100 Myr continuously star- 
forming template with a dust reddening of E(B - V)=0.25) with 
a model where we change the star-formation law to a starburst. 
We find the Schechter parameters to change in the directions that 
are expected for a redder template, as explained above. However 
the differences are insignificant since they are much smaller than 
the statistical errors on the Schechter parameters. 

We stress again that we use a mix of dust amounts in our 
standard analysis to match the UV-continuum slope distribution 
that is measured from the data. Especially for the u-, and g- 
dropouts this puts strong constraints on the combination of the 
age and the amount of dust in the model template, so that we can 
reduce the systematic error to a minimum. 

To justify the assumptions we make regarding the shapes of 
our simulated sources we also estimated the systematic error on 
the LF due to this component. Because we expect similar re- 
sults in the three dropout samples, we only run these simulations 
for the g-dropouts. We inject sources that are 0.05" larger than 
the measured position-dependent PSF, and compare values of 
Moffat parameter B = 3.0 and B = 5.0 with our fiducial B = 4.0 
parameter. We find the following: 

- As expected, a becomes steeper for more extended sources 
(i.e. increasing the Rq or decreasing B), as this causes the 
peak surface brightness to drop. This only significantly af- 
fects the Veff 's at the faint end of the LF. Estimations of or 
range from - 1 .94 < cr < - 1 .36. A similar change is expected 
for the u-, and r- dropout studies. The other Schechter pa- 
rameters, <p* and M*, then also change slightly as they move 
in the direction of the degeneracy in Fig. [9] 

Furthermore we find that, based on the different tem- 
plate spectra, the estimated UV luminosity density varies. The 
ranges are, upon integration down to L < Q.3L* =V in units 
[lO^ergs" 1 Hz 1 Mpc" 3 ], 1.30 < p uv < 1.98 for the u-, 0.92 < 
puv < 1.15 for the g-, and 0.80 < puv < 0.92 for the r-dropouts. 

In order to illustrate the effect that the Eddington bias can 
have on LF estimations we repeat our analysis with slightly 
changed p(m,z). In the completeness simulations we bin the re- 
covered sources by their intrinsic magnitude instead of their re- 
covered magnitude. Doing so leaves the Eddington bias uncor- 
rected. We find that the Schechter fit to the resulting LFs are not 
as good as for our standard analysis, especially for the g-, and r- 
dropouts, where we find that^ nin /dof would be 1.0 and 3.0, re- 
spectively. For the M-dropouts the best-fitting parameters are not 
changed significantly, but for the g-, and r-dropouts the faint-end 
slope steepens while the normalisation of the Schechter function 
decreases. 

Note that we did not account for the presence of Lyman- 
a emission in our simulatio ns, although this lin e clearly con- 
tributes to broadband fluxes. IShaplev et all d2003l ) measured the 
contribution of this line in the spectra of z ~ 3 LBGs and con- 
cluded that only ~ 25% of the sample showed significant Ly 
a emi ssion such that EW(Lyo') > 20A, see also iReddv et al.l 
(HH). Although the relative contribution of this line is thought 
to increase towards higher redshift and fainter continuum lumi- 
nosity, this has not been quantified yet. We will therefore de- 
scribe possible biasses due to the presence of this line in the 
measurement of the LF only qualitatively. 
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If the line appears in emission, it contributes to the flux in 
the CFHT g-band for redshifts of about 2.5 < z < 3.5, in the 
r-band for redshifts of about 3.7 < z < 4.6, in the /-band for 
redshifts of about 4.8 < z < 5.9, and in the z-band for red- 
shifts of about z > 5.8. F ollowing the redshift distributions from 
iHildebrandt et al.l (H009), we expect that the line predominantly 
contributes to the middle band in a two colour selection for the 
m-, and r- dropout samples, moving the source to the upper left in 
Fig- EH This effect causes the effective volumes to rise, thereby 
lowering the LF measurement. If the line indeed gets stronger 
at faint magnitudes, this would bias the LF measurement and 
results in a shallower alpha. In our g-dropout sample the line 
is expected to contribute to the flux in both the r- and /-band, 
depending on the redshift of the particular source. In the case 
where the line falls in the r-band, the effective volumes rise and 
the LF we measured is too high. For the redshifts where the line 
falls in the /-band, the sources would move to the right in Fig.fTOl 
so that the effective volumes are decreased and the LF increases. 
To estimate which effect prevails, the redshift distribution needs 
to be measured spectroscopically. Some objects show Ly-a in 
absorption rather than emission, which would give the opposite 
effect. 

Note that iBouwens et a D (I2007h model the contribution of 
Lya in the measurement of their LF, by using a simple model 
where 33% of the z ~ 4 - 5 LBGs have EW(Ly a )= 50 A, inde- 
pendent of the continuum luminosity. They find that this affects 
the normalisation of the LF by only ~ 10%. However, we stress 
again that the measure of a might be biased if the strength of the 
line depends on the continuum luminosity. 

4.4. The evolving galaxy population 

Although we can measure the UV luminosity density with great 
accuracy, an estimate of the SFR Density depends sensitively 
on the dust ex ti nction correction. Using the prescription from 
IBouwens et al.l d2009l) . we find that the SFR Density shows a 
significant increase by a factor of 4-5 between z ~ 5 and z ~ 3 . 

We find a strong increase in <p* between z ~ 5 and z ~ 3 by 
a factor of 2.5, which is a robust result. The characteristic mag- 
nitude, M* v , however is not very well constrained by our data 
set. This is due to uncertainties in the redshift distributions of 
the source galaxies, as there are no spectroscopic data available. 
Our data are consistent with a non-evolving M^ y ( 20.9) be- 
tween z ~ 5 and z ~ 3. Our data do not show an evolution of the 
faint-end slope, and indicate a 1.6 in this epoch. 

5. Comparison with previous determinations 

Before we compare our results from each of the samples with 
Schechter parameters reported from previous determinations in 
the literature, there are a few things that should be noted. We 
will compare results at redshifts of around 3, 4 and 5. The 
LBGs are generally selected from different surveys and filter 
sets, and therefore intrinsically slightly different galaxies may 
be selected, studied and compared. Also the rest-frame wave- 
length at which the LF is estimated, varies. We compared the 
LFs at 1600A, while e.g. Sa wicki & Tho mpson (2006) measured 
the LF at 1700A, and e.g. lGiavaiisco et al.1 (120041) did their anal- 
ysis at 1500A. Some of the studies described below make use 
of identical or partially overlapping datasets. Our analysis, on 
the contrary, is completely independent from previous determi- 
nations, except for the dust extinction correction in our SFR 
density estimate. Furthermore the errors of several other studies 



could be underestimated because only the Poisson noise com- 
ponent is taken into account, as other noise components (e.g. 
cosmic variance) are difficult to estimate with the typical small 
survey volumes of other surveys. Also note that we compare our 
2-dimensional error ellipses with 1 -dimensional error bars that 
were attained after marginalizing over the two other Schechter 
parameters. This dilutes information on degeneracies between 
Schechter parameters. For these reasons the comparisons in this 
section will be rather qualitative, and are meant to put our re- 
sults into context. Comparisons of the Schechter parameters are 
shown in Table|4] and in Figs.QT|[T2]&[T3] an d will be discussed 
below. 



5.1. Comparison at z=3 

We compare the results from our M-dropout sample with 
several Schechter paramet ers reported in the literature. 
ISawicki & Thompsonl ( |2006|) estimated the LF from the Keck 
deep fields, from t/„G/?-selected star-forming galaxies. Their 
survey area is about 60 times smaller than the CFHTLS-Deep. 
The depth of their observations is slightly deeper than ours, but 
due to our Eddington bias correc tion, we are able t o prob e the 
LFs down to similar magnitudes. iReddv & Steidel (l2009h esti- 
mate a LF at z ~ 3, from 31 spatially independent fields, having 
a total area of about a quarter of ours. Their sampl e contains sev- 
eral th ousands of spectroscopic redshifts at z=2-3. lArnouts et aTl 
(120051) mainly focussed on the galaxy LF at lower redshifts 
(0.2 < z < 1.2) from GALEX-data, but their redshift range was 
extended using 173 galaxies at z ~ 3 from the HDF sample. 
iPoli et al.1 l200lh used the HDF-N, HDF-S and NTT-DF sam- 
ples to estimate the LF in the range 2.5 5= z ^ 3.5. Their sample 
was therefore selected from a very small v olume, which makes 
their results susceptible to cosmic variance. ISteidel et al.l(ll999l) 
pioneered this work, and estimated the UV LF from 0.23 deg 2 
of moderately deep data. Their study was supported by a spec- 
troscopic redshif t samp le. This data is included in the study by 
IReddv & Steidel (120091) . 

Our results agree within the 1 — cr level with the results from 
previous determinations at z ~ 3. Note that the other data points 
lie in the direction of the elongated ellipse, and therefore in the 
direction of the degeneracy we find. 

5.2. Comparison at z=4 

Our g-dropout sample is compar ed with the z ~ 4 sample 
from ISawicki & Thompsonl d2006l) . selected from the GRI fil- 
ter sets on the Keck telescope. Their area is ide ntical to the 
one fr om which the z ~ 3 LF was estimated. ISteidel et al.l 
(fl999l) also estimated the z ~ 4 LF with a similar filter set as 
ISawicki & Thompsonl d2006l) . but did not probe deep enough 
to be able to constrain the faint-end slope a. This parameter 
was set equal to the value found at z = 3, namely a = -1.6. 
lYoshida et al] (120061) presented the LF for 3808 B RV -selected 
LBGs, selected from the Subaru Deep Field project. lOuchi et all 
(120041) selected a z ~ 4 LBG sample from Subaru imaging, 
supp orted by a sample of 8 5 spectroscopically identified ob- 
jects. iGiavalisco et alJ d2004l) used a ~0.09 deg 2 sample from 
the GOODS to estim ate a LF for B45oV606Z850-selected LBGs. 
IBouwens etaD d2007l) used the deep HST ACS fields, includ- 
ing the HUDF and the GOODS, to select a sample of 4671 B- 
dropouts, from which they estimated the U V LF to M1600 ab = 
-16.26. 
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Table 4. Our estimated Schechter parameters compared with values reported in the literature, for the three dropout samples. 



Reference M* uv <p" a 

[lfr 3 Mpc~ 3 ] 

w-dropouts 



This work, Schechter -20.94+JH 4 , L79 -038 -1.65 + ""jj 

This work, Schechter+Magnification -20.84!° : |= 2.11^ : | -1.60+° : ]| 

Reddv & Steidel (2009) -20.97 ± 6.14 1.71 ±0.53 -1.73 ±6.13 



Sawicki & Thompson (2006) 
Arnouts et al. (2005) 

Polietal. (2001) 
Steidel et al. (1999) 


-20 90 +0 - 22 
-21.08 ±0.45 
-20.84 ± 0.37 
-21.04 ±0.15 


1 70 +0 59 

1 ■ ' "-0.32 

1.62 ±0.90 
2.3 
1.4 


-1 43+ 017 

^ -0.09 

-1.47 ±0.21 
-1.37 ±0.19 
-1.60 ±0.13 


g-dropouts 


This work 


-20.84 ± 0.09 


1.36 


+0.23 
-0.20 


-1.56 ± 


0.08 


Bouwens et al. (2007) 


-20.98 ±0.16 


1.3 d 


:0.2 


-1.73 ± 


0.05 


Yoshida et al. (2006) 


-21 14+ 014 

-0J5 


1.46 


+0.41 


-1.82 ± 


0.09 


Sawicki & Thompson (2006) 


-21.0 + °j 


0.85 


-M 

-0.45 


-1 26 +0 - 40 


Giavalisco (2005) 


-21.20 ±6.04 


1.26 i 


:0.03 


-1.64 ± 


0.10 


Ouchi et al. (2004a) 


-21.0 ±0.1 


1.2 i 


:0.2 


-2.2 ± 


0.2 


Steidel et al. (1_999) 


-21.05 


1. 


1 


-1.6 (fi 


xed) 


r-dropouts 


This work 


-20.94+Uf/ 


0.83 


+0.15 
-0.14 


-1 65 +u - (,J 

1,u -0.08 


Bouwens et al. (2007) 


-20.64 ±0.13 


1.0 i 


:0.3 


-1.66 ± 


0.09 


Oesch et al. (2007) 


-20.78 ±0.16 


0.9 d 


:0.3 


-1.54 ± 


0.10 


Iwata et al. (2007) 


-21.28 ±0.38 


0.41 


+0.29 


-1 48 +0 - 38 

i, ^°-0.32 


Yoshida et al. (2006) 


-20 72 +016 


1.23 


-M 

-0.27 


-1.82 (fixed) 


Giavalisco (2005) 


-21.06 ±0.05 


0.83 d 


:0.03 


-1.51 ± 


0.18 


Ouchi et al. (2004a) 


-20.3 ± 0.2 


2.4 d 


: 1.0 


-1.6 (fi 


xed) 



-0.5 



-2.0 



-2.5 



-0.5 



-2.0 



-2.5 



-2.4 
-2.6 
-2.8 
-3.0 
-3.2 
-3.4 
-3.6 
-20.0 
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-21.5 



-22.0 
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Fig. 11. Comparison at z ~ 3. The ellipses represent the 68% 
and 95% confidence contours for different Schechter parame- 
ter combinations, based on our study. The error bars reflect the 
68% confidence limits on the results from previous studies. Our 
results agree within the l-cr level with most of the other deter- 
minations. 



With this data set we have been able to measure the 
Schechter parameters for the g-dropout sample with very high 
statistical accuracy. Note however that several systematic uncer- 
tainties, which we comment upon in Sect. 14.31 are not included 
in our error ellipses. Our results agree within the l-cr level 
with many of the z ~ 4 results in the literature. However, there is 
still some tension in measurements of the faint-end slope a. The 
characteristic magnitude, M^ v , we measure, is slightly fainter 
than the values we found in the literature. 



Fig. 12. Similar as Fig. [TT] but now a comparison at z ~ 4. 



5.3. Comparison at z=5 

At z ~ 5, llwata et al] (120071) reported the UV LF from a combi- 
nation of H DF and Subaru imag es, totalling a survey area about 
1/9 of ours. lOesch et"ai] (l2007h based their study on approxi- 
mat ely 100 LBGs from very deep ACS and NICMOS imag- 
ing. lYoshida et al.l (120061) also defined a z ~ 5 sample from 
their observations , combi n ing Vi'zf and Ri ' z' sele cted objects, 
as did lOuchi etaHd2004al) . lGiavaiisco et al.ld2004l) selected 275 
y606'77 5Z850 LBGs to estimate a z ~ 5 UV LF. iBouwens et al.l 
(2007) also measured a sample of 1416 V-dropouts from their 
deep HST ACS sample, which resulted in an estimation of the 
UV LF down to M l600AB = -17.16. 

Similar to the Schechter parameters found for z ~ 4, there is 
a large discrepancy in the literature for the Schechter parameters 
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Fig. 13. Similar as Fig. [TT] but now a comparison at z ~ 5. 
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Fig. 14. psfr as a function of redshift and cosmic time. Dotted: 
SFRD uncorrected for dust. Solid: the dust corrected SFRD, 
where we use a lumino sity dependent dust correction factor from 
iBouwens et al.l d2009l) . Note that we included both random and 
systematic errors in our dust correction, unlike most other stud- 
ies. 



at z ~ 5. The statistical uncertainties in the Schechter parameters 
is very small for our r-dropout sample. Note however that several 
systematic uncertainties are not included in these error ellipses, 
see Sect. 14. 31 Our results agree reasonably well, within the 1 ~ cr 
level, with many previous determinations at z ~ 5. 



5.4. Comparison of the SFR density 

In Fig. [I4]we com pare the SFR density values given in Table[3]to 
values reported bv lSchiminovich et al.l (120051) . who made use of 
low -z GALEX data iReddv & Steidef (120091) at intermediate z, 
and lBouwens et all d2009l) at high z. The uncorrected SFRDs are 
in good agreement with each other and show a smooth redshift 
evolution. However, it is clear that the dust correction is the ma- 
jor uncertainty because of the age-dust degen eracy. We use the 
same dust correction as IBouwens et al. I J2009I) and also include 
systematic uncertainties in the error bars. 



6. Summary & Conclusions 

In this paper we use the CFHT Legacy Survey Deep fields to 
estimate the UV Luminosity Functions of the largest u, g-, and 
r-dropouts samples to date. As our samples are all extracted from 
the same dataset this study is ideally suited to study a time evo- 
lution of the luminosity function in the redshift regime z = 3-5. 
Thanks to the large volumes we probe with our 4 square de- 
gree survey, cosmic variance plays a negligible role in our anal- 
ysis. We are now able to study the bright end of the luminos- 
ity function with unprecedented accuracy. Furthermore, given 
the depth of the stacked MegaCam images, we probe the faint 
end of the luminosity function with comparable precision as the 
deepest ground based surveys have done before. This unique 
combination gives us the opportunity to not only estimate the 
Schechter parameters for the different luminosity functions, but 
also to study a possible deviation from this commonly used fit. 

In the m-, and g-dropout samples of our survey we are able 
to measure the UV continuum slope directly from the data. This 
allows us to simulate sources that have the same distribution of 
UV slopes, which is important for an accurate estimate of the 
Schechter parameters. 

We find the faint-end slope a to not evolve significantly in 
the redshift range we probe, and to have a value of around -1.6. 
This parameter however, is not very strongly constrained by our 
ground based survey, as this parameter depends on some of the 
assumptions made. 

We do not find a significant evolution in MJj v , and argue 
that this might be due to insufficient knowledge of the redshift 
distribution of the source galaxies. The conversion of apparent 
to absolute magnitudes depends strongly on these distributions, 
and an uncertainty in the distance modulus directly propagates 
into an equal uncertainty in M^ v . This parameter is therefore 
poorly constrained by this study, until a more reliable redshift 
distribution is available. 

We find a strong evolution in <p*, which we argue to be sig- 
nificant. The normalisation of the LBG density, <p*, increases by 
a factor of ~ 2.5 from z ~ 5 to z ~ 3, an increase that cannot 
be explained by any change in the assumptions tested. We there- 
fore conclude that the UV luminosity density is increasing in the 
corresponding epoch, in a way that does not strongly differ with 
magnitude. 

The SFR Density does increase significantly, by a factor of ~ 
3, between z ~ 5 and z ~ 4. We find a smaller, but less significant 
increase between z ~ 4 and z ~ 3. 

With our 4 square degree survey we probe densities that are 
at least four times lower than any of the studies we compared 
our results to. We find a substantial deviation from the Schechter 
function at the bright end for the M-dropouts, where the LBG 
densities are very low. We find that the deviation can be at- 
tributed to magnification effects that arise from inhomogeneities 
in the matter distribution between the LBGs and the observer. 
We fit an improved Schechter function that is corrected for mag- 
nification and find that the quality of the fit improves signifi- 
cantly. Intrinsically the distribution of luminosities does there- 
fore not deviate significantly from a Schechter model. With this 
data set we have been the first to be able to measure a hint of this 
magnification imprint on a z ~ 3 LBG sample. 
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Fig. 10. The grey points represent the colours of 10,000 objects in field 1 of the CFHTLS Deep. The black boxes are defined by 
Eqs. Q][3]and are used to select M-(top panels), g-(middle panels), and r-(bottom panels) dropouts. The coloured tracks represent 
the colours of a template galaxy as a function of z. They are evaluated at intervals of Az=0. 1 . The symbols mark redshifts from 
2.5 to 3.5, from 3.2 to 4.2 and from 4.2 to 5.2 for u-, g-, and r-dropouts respectively. The green curve with x-symbols represents 
a reference model, a 100 Myr old continuously star-forming template with a dust attenuation of E(B - V)=0.25. In the left panels 
we consider redder (bluer) templates by increasing (decreasing) the template age, see the red (blue) curves and the A (O)-symbols. 
In the right panels we consider redder (bluer) templates by increasing (decreasing) the amount of dust in the template, see the red 
(blue) curves and the A (O)-symbols. Note that we use a distribution of dust in our fiducial analysis (see Sect. 13. lb . 



